Initial state dependence in multi-electron threshold ionization of atoms 
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It is shown that the geometry of multi-electron threshold ionization in atoms depends on the initial 
configuration of bound electrons. The reason for this behavior is found in the stability properties of 
the classical fixed point of the equations of motion for multiple threshold fragmentation. Specifically 
for three-electron break-up, apart from the symmetric triangular configuration also a break-up 
of lower symmetry in the form of a T-shape can occur, as we demonstrate by calculating triple 
photoionization for the lithium ground and first excited states. We predict the electron break-up 
geometry for threshold fragmentation experiments. 
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Three-body Coulomb dynamics, in particular two- 
electron atoms, are very well studied in the energy regime 
of single as well as double ionization [l], [3j S 01 ■ Much 
less is known about correlated dynamics in four body 
Coulomb systems, more precisely on differential observ- 
ables for fragmentation of a three-electron atom in its 
nucleus and all electrons 0, @, 0] ■ A recent experiment 
provides for the first time detailed information in terms 
of differential cross sections on the angular and energetic 
break-up parameters of three electrons following impact 
double ionization of Helium Q . For small excess energies 
E (each continuum electron carries away about 9 eV en- 
ergy) , it was found that the electrons form an equilateral 
triangle upon breaking away from the nucleus. This is 
expected in accordance with Wannier's theory Q, quan- 
tified for three electrons in There, it is shown that 
the fixed point of classical dynamics, through which full 
fragmentation near threshold E = Q should proceed, is 
given for a three-electron atom by an equilateral triangle 
with the nucleus in the center and the electrons at the 



In two-clcctron atoms the corresponding fixed point 
implies a collinear escape of the electrons in opposite di- 
rections [^, [lH . The normal mode vibration about this 
collinear configuration is stable, which is in marked con- 
trast to the three-electron case, where the triangular con- 
figuration is linked to two unstable, degenerate normal 



modes 
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We will show that as a consequence of this property the 
preferred final geometry of the three escaping electrons 
becomes initial state dependent and can change between 
an equilateral triangle and a less symmetric T-shaped 
escape. While the former is realized, e.g., in electron im- 
pact double ionization of helium, the latter should be seen 



in triple photoionization of lithium. These are only two 
prominent examples, the general pattern and the reason 
for it will be detailed below. 

Due to the scaling of the Coulomb potential, states of 
finite total angular momentum L will all behave like the 
L = state close to threshold which is therefore sufficient 
to consider [isj . In hyperspherical coordinates with the 
radial variable w instead of the hyperradius r = w'^, the 
Hamiltonian for a three-electron atom with total angular 
momentum L = reads 



8w^ 



(1) 



where $7 = (ai, a2, 9i, ^2)^ contains all angular variables 
describing the positions of the electrons on the hypcr- 
sphcrc of radius r and A, the so called grand angular 
momentum operator , is a function of Jl and all conju- 
gate momenta. The total Coulomb interaction V = C /r 
acquires in this form simply an angular dependent charge 
C{n). In terms of the familiar vectors ri,r2,r3, point- 
ing from the nucleus to each electron, the hyperspherical 
coordinates are given by 



r 




(2a) 


ai 


= arctan(ri/r2) 


(2b) 


a2 


= arctan(r2e/r3) 


(2c) 


01 


= arctan[ri • r2/(rir2)] 


(2d) 


02 


= arctan[ri • r3/(rir3)] , 


(2e) 
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where r2e = [rf + t"!)^^^- 

Threshold dynamics is governed by motion along the 
normal modes about the fixed point of the Hamiltonian 
H — {h — E)w'^ The special form of H ensures that 
the dynamics remains regular approaching the fixed point 
radially, i.e., w — > w* = 0, while Q,* at the fixed point 
is defined through VnC(Jl)|f2=o' = 0. The equations 
of motion can be expressed as a system of first order 
differential equations (ODE) T = QVrH for the phase 
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space vector T = (p^, P^, w, fJ)''', with 

a block matrix composed from and unity matrices of 
dimension f x f [16[, where / is the number of degrees 
of freedom, here / = 5. Since the differential equations 
are still singular at the fixed point {w* , ft*), a change of 
the momentum variables Pq conjugate to Q is needed 

Pu^^^PuJw (4) 

as well as a new time variable r related to the origi- 
nal time t conjugate to the Hamiltonian Eq. ^ through 
dt = w^dr. Finally, the normal modes can be obtained 
from the modified ODE, dj/dr = QV^H by diagonal- 
izing the matrix d'^QH /{d'ydj)\-^=-^-' , where 7 refers to 
the new phase space variables with the (non-canonical) 
momenta from Eq. (|4|). The eigenvalues are the Lia- 
punov exponents Xj and in the normal mode basis {uj} 
threshold dynamics assumes an oscillator-like form of 
duj{T) = cxp{XjT)Suj(0) with the unit vectors iij 

iij = Suj{0)/\Suj{0)\ (5) 

defining the normal mode basis. The <57(t) =7 — 7* are 
excursions of 7 from their fixed point values 7* and are 
expressed as a linear combination of the 6uj (r) . 

We recall briefly the familiar three-body break-up in 
a two-electron atom with hyperradius r2e and the angles 
ai,9i defined as in Eq. ([2]). The charge corresponding to 
C(f2) in Eq. ([2]) is for the two-electron problem 



>-^2eV"l7 '^IJ — : r ^ -T ■ 

smai cosai (f - sin(2ai) cos 6*1) 2 

(6) 

The fixed point analysis reveals a pair of unstable A1/2 = 
— Ao ± A and stable A3/4 = —ujq i iuj Liapunov expo- 
nents (Ao,a;o,A,w > 0) with a shift (— Aq and —luq, re- 
spectively), compared to standard symplcctic dynamics. 
The shift formally arises through the noncanonical trans- 
formation of the momentum variables Eq. (|4|) necessary 
to obtain normal mode motion about the singular fixed 
point w* = 0. 

The resulting eigenvectors iii reveal orthogonal motion 
along 9i and ai, i.e., any phase space vector 157,^^ (r) de- 
scribing linearized motion in the subspace spanned by 
, Q^i, can be expressed as a linear combination of two 
eigenvectors ^70,^ = a-y^ exp(Air)'Ui -I- b-y^ exp(A2T)'U2. 
An analogous relation holds for linearized motion in the 
subspace , 0i , realized through two different eigenvec- 
tors, (57g^ = a-yg exp(A3T)il3 -I- b^g exp(A4r)'U4. Hence, 
^leii''') ■ ^laii''') = at all times r. The coincidence of 
the eigenspaces with the respective dynamics of 9i and ai 
has the important consequence that the fixed point value 
01 = TT is preserved through its relation to the stable 



TABLE I: The two eigenvectors Ui belonging to the positive 
Liapunov exponent A > in the basis of the phase space 
variables from Eq. ((2} 



basis 


Ui 


U3 




-8.22 X 10"^ 


2.67 X 10"^ 




-2.97 X 10"^ 


9.64 X 10"^ 




-1.81 X 10~^ 







-9.83 X 10"^ 





661 


7.12 X 10"^ 







6.57 X 10"^ 


-4.46 X 10"^ 


502 


3.33 X 10"^ 


7.40 X 10-3 


Spe2 


-2.75 X 10"* 


8.92 X 10-3 



eigenmode while all energy sharings can occur through 
the relation of the unstable eigenmode with ai. 

This coincidence of normal modes with subspaces of 
obscrvablcs is special to two-clcctron atoms and does not 
hold for more electrons. Moreover, for three-electron 
atoms, a new feature emerges in threshold dynamics, 
namely the existence of two degenerate pairs of unstable 
normal modes with Liapunov exponents X1/2 = A3/4 ~ 
— Ao ± A. Note that the fragmentation dynamics close 
to the fixed point will take place in the phase space of 
the two unstable normal modes with an arbitrary lin- 
ear combination of the two eigenvectors ui and U3 be- 
longing to the two (equal) positive Liapunov exponents 
Ai = A3 = A+, 

Sj{t) = exp(A+T)(ciiti -I- C3-U3) , (7) 

which allows for flexibility in the four-body break-up, as 
we will see. 

One sees directly from Table I that normal mode dy- 
namics about the fixed point for geometrical angles 9i and 
hyperanglcs ai is not separated as for the two electron 
case, and even more importantly, all phase space vari- 
ables are linked to the unstable normal modes. Hence, 
the fixed point geometry docs not provide necessarily a 
preference for the final angles of the electrons. On the 
other hand, this opens the way for the initial state to 
have an influence on the final observables, even close to 
threshold. 

Relevant for threshold ionization is the spatial electron 
distribution at the time (we label it r = 0) when all elec- 
trons to be ionized have received enough energy (through 
collisions) to leave the atom. We call this distribution the 
transient threshold configuration (TTC). 

In a two electron-atom, the necessary energy transfer 
between the two electrons leading to three-body fragmen- 
tation happens through a single collision. At this time 
T = both electrons are naturally close together, so that 
Sai{0) w holds and defines the TTC, independent of 
the initial bound electron configuration. 
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In a three-electron atom, the situation is more com- 
plicated since at least two collisions are necessary to dis- 
tribute the energy among the electrons so that all of them 
can escape. For triple photo-ionization of Lithium we 
know from classical calculations that the Is photo elec- 
tron (3), which has absorbed the photon initially, collides 
immediately with the other Is electron (2) and subse- 
quently (about 60 attoseconds later) either electron 2 or 
the photo electron itself collides with the 2s electron (1). 
This can be expressed with the two collision sequences 
si = (32,21) and S2 = (32,31) [2lj. The time delay of 
the second collision, respectively, is due to the "distance" 
of the 2s shell from the Is shell. It leaves an asymmet- 
ric situation after the second collision when the transient 
threshold configuration is reached at t = 0. While the 
two electrons participating in the last collision are close 
to each other, the third one is further away. Concen- 
trating on the collision sequence (32,21), this implies in 
terms of distances to the nucleus 



n « r2 7^ rs . 



(8) 



Hence, Sai{0) w while Sa2{0) ^ 0. This can be easily 
accommodated with a suitable linear combination choos- 
ing, e.g., ci ~ 1/a and C2 = —1/6 in Eq. ([71), where a 
and h are the coefficients of 8a\ with respect to u\ and 
M3 (see first row of Table I) . 

The scenario described is still within the overall pic- 
ture of threshold break-up in the spirit of Wannier and 
therefore in accordance with the power law a oc [EjEt^)'^ 
for the dependence of the total break-up cross section 
on the excess energy E [13, 18, [l^. The exponent 
a = (Ai A3)/(2A^) = 2.162, where A^ = 2.5088 is the 
radial Liapunov exponent and Ai = A3 = 5.4240. Yet, 
the TTC with unequal distances of the electrons to the 
nucleus breaks the complete symmetry among the elec- 
trons. This in turn, leads to a preferred final geometry of 
the three electrons since essentially with the third elec- 
tron being far away, the geometry determination reduces 
to that of two electrons escaping in an equivalent way as 
in the two-electron break-up, i.e., back-to-back with an 
angle oi d\ ~ tt. The question remains if there is a pre- 
ferred angle between this electron pair escaping along a 
line and the third electron. Assuming for simplicity that 
rs ^ ri K r2 and therefore sin Q!2 ~ holds (the oppo- 
site case would lead to the same result), wc can expand 
C{Q) in Eq. ([T]) in powers of rj = a2 
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FIG. 1: Probability to find the angle 9 between two electrons 
for triple ionization from the initial Li(ls2s^) state for excess 
energies above the triple ionization threshold as indicated. 



with 



Co = 




(10a) 


Cl = 


2-Z 


(10b) 


C2 = 


cos 02 sin ai + cos 612 cos ai 


(10c) 


C3 = 


|(cos6'2 sinai)^ -|- |(cos6'i2 


cos Q!i)^ 




+(l-2Z)/6, 


(lOd) 


S12 = 


02 — Oi. To lowest order in 


77, the problem 



c„r] ' 



(9) 



n=0 



to find a stable configuration is that of the two-electron 
system (here with nuclear charge Z = 3) with the well 
known solution 0** = tt and al* = 7r/4 [§]. These val- 
ues minimize C2 for any value 02- Its value 02* = 7r/2 is 
determined from dcs/d02 ~ which is a stable solution. 
This constitutes the T-shape of the three electrons as a 
preferred asymptotically (77 —>■ 0) stable geometrical con- 
figuration within the globally unstable two-dimensional 
subspace spanned by ui^u^. The fragmentation dynam- 
ics in the sub space is completely degenerate - hence, we 
have looked on a "higher order" correction, which would 
give a preference within the degenerate subspace: that is 
the asymptotically stable T-shape. Indeed, this configu- 
ration was found numerically close to threshold in triple 
photo ionization [20| . 

One can double check this new insight into the role of 
the TTC of the electrons in the presence of degenerate 
Liapunov exponents by considering initial configurations 
which lead to different transient configurations. 

For an initial excited state Li(ls2s^) the two collisions 
in which energy is transferred from the (Is) photo elec- 
tron to the two 2s electrons happen close in time (and 
therefore in space). Hence, the transient threshold config- 
uration after the second collision at r = is ri « r2 ~ r3 
close to the fixed point and therefore giving preference 
to the symmetric break-up with a dominant angle 120° 
between the electrons, similar to electron impact double 
ionization of Helium. 

Numerical results confirm this prediction as seen from 
Fig. [1] where the probability to find an angle between 
two electrons in triple photo ionization of the initial 
Li(ls2s^) is shown for different excess energies. Clearly, 
the most likely angle is cos 6* — —1/2 corresponding to 
120°, indicated by the thin vertical line. 
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FIG. 2: Same as Fig. [TJ but for triple ionization from the 
initial Li(ls'^2s) ground state for excess energies above the 
triple ionization threshold as indicated. 
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FIG. 3; Same as Fig. [TJ but with a 2s electron as photo elec- 
tron. 



On the other hand, for the ground state as initial state 
the TTC is asymmetric as previously described and we 
expect a final T-shape geometry with peaks at 90° and 
180°. This is indeed the case, as can be seen in Fig. [21 
As a final test we may use again the excited initial state 
Li(ls2s^) but take the 2s electron as the photo electron. 
This process is due to the smaller dipole coupling by 
more than an order of magnitude suppressed compared to 
the ionization with the Is electron absorbing the photon. 
However, here we use this only as an illustration for the 
initial state dependence of threshold ionization. Accord- 
ing to our reasoning we have in this case (although it is 
the same initial state as before) a different TTC since the 
first collision of the photo electron 1 happens with the Is 
electron 2 while later on the collision of the Is electron 
2 with the 2s electron 3 will take place. Consequently, 
for the TTC fa r2 ^ ri holds which is structurally 
identical to Eq. ([S]), the situation when a final T-shape 
geometry appears. As can be seen in Fig. [3] the T-shape 
indeed also emerges if the photo electron comes from the 
2s shell. 

To summarize, we predict that for three electron break- 
up near threshold two preferred geometrical patterns for 
the electrons exist, namely an equilateral triangle and 
a T-shape. Which of them is realized depends on the 
transient threshold configuration, that is the spatial dis- 
tribution of the electrons at the time when the energy 
among them is distributed such that all of them can es- 
cape. The transient threshold configuration is strongly 
influenced by the initial state of the electrons, as has 
been discussed in detail for the case of lithium. A com- 
plete overview of the three-electron break-up pattern for 



TABLE II: The relation of excitation process, initial configu- 
ration and geometry of the electrons for various three-electron 
break-up processes close to threshold. 



impact 


initial 


break-up 


by 


configuration 


geometry 


electron 


Is^ (He) 


A 


electron 


ls2s (He) 


_L 


photon 


ls^2s (Li) 


_L 


photon 


Is2s2 (Li) 


A 


photon 


Is2s3s (Li) 


_L 



three electron systems is provided in Table [TTl When- 
ever the two electrons to which energy is transferred in 
the course of the fragmentation of the atom are in the 
same shell, a symmetric triangular geometry is expected. 
If the two electrons are from different shells, we expect 
a T-shape. Recent experimental results on electron im- 
pact double ionization of helium Q are consistent with 
our predictions but do not provide sufficient information 
to consider our prediction as experimentally already con- 
firmed. Based on the present consideration it is also pos- 
sible to predict the break-up geometries and their real- 
ization dependent on the initial state for more than three 
electrons. 
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